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I. INTRODUCTION 

Quantum chrome-dynamics (QCD) at non-zero tem- 
perature is an intensely studied field. At very high tem- 
peratures the coupling constant is small and perturbation 
theory can be developed. However, due to the chromo- 
magnetic sector of the theory, perturbation theory ex- 
plodes already at a few-loop level 0, 0, and is hence 
only applicable at academically high temperatures Q. 
The region of intermediate temperatures is of much big- 
ger interest. Both the restoration of chiral symmetry 
and the deconfincment are believed to take place in this 
region. In the presence of fermions it is still unclear if 
there is a confinement-deconfinement phase transition or 
just a smooth crossover between the two phases. In any 
case QCD is in the deconfined plasma phase at very high 
temperatures. 

At finite temperature gluons obey periodic and 
fermions obey anti-periodic boundary conditions in imag- 
inary time. This property leads to the quantized Mat- 
subara frequencies. They are even multiples of ttT for 
gluons, i.e. ujk = 2kitT, and odd multiples of 7rT for 
quarks, i.e. = (2k + l)7rT. So while gluons have a 
zero mode this is not the case for quarks. This fact has 
direct and important influence on the IR behavior of the 
two contributions to the effective action. 

At the tree level very heavy modes decouple from a 
theory at high temperatures. This is called dimensional 
reduction |5[ since the heavy modes are simultaneously 
the time-dependent ones. Neglecting all modes except 
the zero Matsubara frequencies leaves a 3D static theory 

- + ^ Wt - - n + m b A\f] , (i) 

which only contains the static gluonic modes with the 
coupling constant g 2 3 = g 2 T . Since the energy can never 
vanish for fermions they decouple completely. 

The long-range forces mediated by the static gluons 
lead to the IR divergencies, because in strict perturbation 
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theory they are massless. Fermions do not cause any 
IR problems, since they do not have zero modes even if 
they are massless, which is the case we consider here. 
Nevertheless, fermions change the effective action in a 
drastic way as compared to the pure glue case, since their 
presence changes the symmetry of the action with respect 
to the center-of-group gauge transformations. 

The tree-level action (TjQl has certainly insufficient accu- 
racy to study field fluctuations at high but not infinitely 
high temperatures. Once one includes quantum correc- 
tions both the fermions and all the non-zero Matsubara 
modes of the gluons show up in the loops. 

Effective theories resulting from quantum corrections 
for different energy scales, T, gT and g 2 T have been con- 
structed in 0, Q . The parameters in the effective theo- 
ries are obtained by matching the correlation functions 
between the effective and the actual theory as functions 
of the parameters of the original theory. The question of 
color-conductivity and transport properties of the plasma 
has been addressed in Q. A heat kernel approach for 
Yang-Mills theories has been used in [j| , a constraint ef- 
fective potential for the Polyakov loop has been studied 
in [lfj and spatial variations of the Polyakov loop have 
been investigated in [Tl|. 

In the pure Yang-Mills theory the center symme- 
try plays a crucial part in the description of the the 
confinement-deconfinement phase transition Q, IT^L Il3l 
fl4| . The latter is usually characterized by an order pa- 
rameter which is the average of the trace of the so-called 
Polyakov line: 

/ ri/T \ 
P(x) = P exp i / dx4 A4 . 

The order parameter < TrP > is zero in the confined 
phase below the critical temperature and assumes a non- 
zero value in the deconfined phase above the critical 
temperature. The Polyakov line is not invariant under 
gauge transformations belonging to the gauge group cen- 
ter. One hence concludes that if < TrP >= then the 
Z(N C ) symmetry is manifest. This situation describes 
confinement. If for any reason < TrP >^ then the 
symmetry must have been broken. This corresponds to 
the deconfmedphase. 

The 1-loop [1 E E3 and 2-loop [T^ potential ener- 
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gies as functions of A4 are known. They are periodic 
functions of the eigenvalues of A4 in the adjoint repre- 
sentation with period 2ttT. This reflects the symmetry 
of the Z(N C ) vacua. The curvature of the potential at its 
minima gives the leading order Debye mass for 'electric' 
gluons. The zero energy minima of the potential corre- 
spond to quantized values of A4 or center group values for 
the Polyakov line, where TrP 7^ 0. At high temperatures 
the system oscillates around one of these minima. At 
lower temperatures, however, the fluctuations around the 
minimum increase and eventually the system undergoes 
a phase transition to < TrP >= 0. At the same time, 
one expects that near the phase transition point the fluc- 
tuations are long-range. To study those fluctuations, one 
needs an effective low-momenta theory which, however, 
does not assume that the A4 component is small. 

Let us formulate the problem more mathematically. 
Nonzero temperatures explicitly break the 4D Euclidean 
symmetry of the theory down to the 3-D Euclidean sym- 
metry so that the spatial Ai and time A4 components 
of the Yang-Mills field play different roles and should 
be treated differently. One can always choose a gauge 
where A4 is time- independent. Taking Aj(a;) to be 
static is not a restriction of any kind on the fields but 
merely a convenient gauge choice, and we shall imply 
this gauge throughout the paper. [It is also a possible 
gauge choice at T = but in that limiting case it is 
unnatural as one usually wishes to preserve the AD sym- 
metry] As to the spatial components Ai(x,t), they are, 
generally speaking, time-dependent, although periodic in 
the time direction. Putting the components Ai to zero 
is a gauge non-invariant restriction on the fields since 
any time-independent gauge transformation will gener- 
ate a nonzero Aj. Therefore, the spatial derivatives of the 
Polyakov line in the gauge-invariant effective action can 
only appear as covariant derivatives including a nonzero 
Ai field. 

In ^5 we calculated the 1-loop kinetic energy for the 
eigenvalues of the Polyakov line, integrating over gluon 
and ghost fluctuations. See also for a summary. In 
this work we are interested in obtaining the effect of 
quarks on that kinetic energy as well. We use a back- 
ground field method for the gluons and evaluate the 1- 
loop action through a functional determinant formalism. 
In particular we assume the background fields to vary 
slowly but the A4 component is allowed to have an arbi- 
trary amplitude. We integrate out fast varying quantum 
fluctuations about them by making an expansion in spa- 
tial covariant derivatives. This method was originally 
developed in 20] for zero temperature QCD. 

This corresponds to summing up all powers of A4 but 
where their momenta are restricted to p < T reflecting 
the long-range behavior of the plasma phase. As we said, 
we choose a static gauge for A^x). This gauge choice 
does not prevent Ai(x) from being time dependent. Since 
Ai(x4,x) is periodic in time, its time derivative is given 
by the Matsubara frequencies cuk = 2irkT being 0{T) 
for any fc ^ 0. Since we are interested in low momenta 



fluctuations, p < T, it is consistent to restrict oneself to 
the zero Matsubara frequency of the background field, 
i.e. to the static Ai(x). 

We expect that our results are suitable to study the 
correlation functions of the Polyakov line not too far 
from the transition point where it experiences fluctua- 
tions that are large in amplitude but presumably mainly 
long ranged. The results may be of some help for study- 
ing quantum weights of semiclassical objects, such as 
dyons (HH,!^) or calorons [23| . 

The effective action contains a contribution from the 
gluons and from the fermions. The former part, namely 
the pure Yang-Mills effective action, was obtained by the 
authors in Although, as discussed above, there is no 
center symmetry for the fermions, it is still instructive to 
see their effect on the effective action for the Polyakov 
line. In particular we work with the gauge group SU(2), 
We consider a general electric field but restrict ourselves 
to a magnetic field parallel to A4, £?J. The expected 
result for the effective action, which is the sum of the 
tree-level and 1-loop actions, is hence: 



=J ^ [-t 3 v f (aD 



-E?F?\Al) 



■{B}fH^{Al) 



(2) 



with the electric and magnetic fields 

E? = DfA\ - A a t = d l A a 4 + e acb A c t A\ - A?, (3) 



m = ^e ijk (d j At-d k A] + e abc A b j At). (4) 



The first term is the potential energy and the remaining 
terms are the kinetic energy contributions in the color- 
electric and color-magnetic sector. The objective of this 
paper is to find these functions. 



II. THE QCD ACTION AT FINITE 
TEMPERATURE 

The basics about Yang-Mills theory at finite temper- 
ature were discussed in [Ig. The (Yang-Mills) action of 
gluons at finite temperature is given by 



d,X4 



(3=^, (5) 



where the gluon fields obey periodic boundary conditions 
in the temporal direction, i.e. 



4^(0, x) = Af,{P,x). 



(6) 



Because of the compactified time direction there is a 
group of special gauge transformations which transform 
the gluon fields in the usual way as 



A„ 



UA^U- 



(7) 
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and which preserve the periodicity condition JfjJ), but 
which are not periodic themselves: 

U(0,x)=z k U{p,x). (8) 

Here z k is an element of the center group Z(N C ): 

Zk ^ e 2mk/N c fcg{0,iV c -l}. (9) 

The Yang-Mills action is invariant under this gauge 
transformation, but the Polyakov line is not. It trans- 
forms as 



( f 1/T 
P(x) = P exp I i J dtA 4 



z^P(x). (10) 



From this property one sees immediately that the Z(N C ) 
symmetry implies < TrP >= while it must be (spon- 
taneously) broken if < TrP >^ 0. 

Including Nf quarks with mass mj the full QCD action 
at finite temperature becomes 



S = dxi 
Jo 



d 3 3 



(11) 



N f 



V>} + im f) ">Pf 



where the Dirac operator is given by 



(12) 



Here the T a are the generators of SU (N) in their funda- 
mental representation, they are half the Pauli matrices 
for SU(2). The fermions in eq. obey anti-periodic 
boundary conditions 



ipf(0,x) = -ipf(j3,x). 



(13) 



This property is, however, not preserved by the Z(N C ) 
gauge transformation (jSJ. Specifically the fermions 
transform as 



(14) 



^(0, a:) = U(0,x)i> f (0,x) = z k U{(3,x)^ f (0,x) 
= -ZkU(j3,x)<ipf{j3,x) =-ZkipJ(0,x). 

Hence in the presence of fermions the Z(N C ) symmetry 
gets explicitly broken and the Polyakov line ceases to 
serve as an exact order parameter for the theory, since 
< TrP >^ for all temperatures. Nevertheless, even in 
the presence of massless fermions the Polyakov line might 
provide useful information near the critical temperature 



III. ONE LOOP QUANTUM ACTION 

The partition function of QCD in its Euclidean invari- 
ant form is given by 



Nf 

Z(A,ip f ,ip f ) = ]T / DA Dip f Dtfj^ exp / d^x (15) 
f=1 J 



N f 



We use the background field method for the gluon fields, 
where we decompose them into background fields and 
quantum fluctuations around them which we assume to 
be small: 



An — Afi + a^. 



(16) 



In this work we are interested in a 1-loop effective theory 
for the background A fields. This corresponds to an ex- 
pansion of the action around the background gluon fields 
to quadratic order in the quantum fluctuations a„. The 
one loop expansion of the gluon Lagrangian is: 



1 9 , 1 



—Ff,JA) = ^ — rFf,JA) 



V(M) 

-^±- )D M)F^{A) ai 

a a W ab a b + 

2g 2 (M) " u 



(17) 



where 



W, 



and 



ab . 



[D\A)] ab 5^[DM ab -2r b F^(A), (18) 



Df{A) 



d»6 



ab 



£acb j^c 
J l 



(19) 



is the covariant derivative in the background field in the 
adjoint representation. The second term in eq. i|17|l . 
which is linear in a„, is zero if the background field obeys 
the equation of motion. In the fcrmionic Lagrangian the 
quarks couple to the gluon fields in the usual minimal, 
i.e. linear way. Hence the expansion is just 



(20) 



where X 7 (A) is the covariant derivative of the background 
field in the fundamental representation. The second term 
in eq. I|20(l contributes at the 2-loop level which we do not 
consider here. 

The quadratic form W® b in eq. l(T%|) is degenerate: it 
has an infinite number of zero modes which are the in- 
finitesimal gauge transformations = Z3° b A h . In order 
to remove this degeneracy one has to fix the gauge for 
these fluctuations. We choose the background Lorenz 
gauge D ll (A)a fl = [2^. The second term in eq. l(T%|) 
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cancels out but the Faddeev-Popov ghost determinant 
arises which again can be expressed as a Grassmann in- 
tegral over ghost fields. 

The 1-loop partition function thus becomes 



Z(A) 



Da D\ D\ + Dip / Dip / exp i / d 4 x 



(21) 



In this paper we study the contribution of the fermion 
functional determinant to that effective action. We use 
again the technique of the covariant derivative expansion. 
In addition, we will work in the chiral limit throughout, 
i.e. we set rrif — 0. The main difference to the gluon cal- 
culation is that in the case of fermions we are dealing with 



1 



2g*(M) 

where x, x + are ghost fields 
i Tp j = iyi + im f = i jz 



N f 

K <l < - x +a D l x a ) +Y^ipfip i>i 

/=1 



T a A a n 



(22) 



is the massive Dirac operator in the fundamental repre- 
sentation, and 



1 



5(mJ J d^ u {A)F^(A) 



(23) 



is the action of the background gluon fields. 

Integrating out the quarks, ghosts and the quantum 
fluctuations of the gluons leaves us with the desired ef- 
fective theory for the background A fields: 



z(A) 



(detW) 



-1/2 



det(-L> 2 ) JJdet(i^j) 



-1/2 



logdet (-D 2 



(24) 
(25) 



so that the 1-loop action is 
Si-loop = log (detW)' 

N f 

+ ^logdet (ip f ) . 
/=i 

Since the A are the only gluon fields left we will omit 
the bar from now on. So far the background field has 
been kept arbitrary. One has, however, the gauge free- 
dom to choose the A±(x) fields to be static. The spatial 
gluon components are generally time dependent. Since 
Ai(x4,x) is periodic in time, its time derivative is given 
by the Matsubara frequencies ojk = 2nkT being 0(T) for 
any k ^ 0. Since we are interested in low momenta fluc- 
tuations, p < T, we shall restrict ourselves to the zero 
Matsubara frequency of the background field, i.e. to the 
static Ai(x). 

The operators in the ghost and gluon functional de- 
terminants, D 2 and W, are matrices in the adjoint rep- 
resentation of the color group, and they are built from 
covariant derivatives and the field strength only. We used 
this fact in 01 to make an expansion of the 1-loop pure 
Yang-Mills action in powers of Di. Since the (static) 
electric field is given by Ef = Df b A\ and the magnetic 
field by B% = leyfcJVJ = \e ijk e cad [D t ,D 3 ] cd we obtained 
an effective action for the background A4 fields in terms 
of electric and magnetic fields. 



operators in the fundamental representation and that we 
do not expect Z(N C ) symmetric results. 



IV. THE FERMIONIC FUNCTIONAL 
DETERMINANT 

Throughout this paper we will be working with Eu- 
clidean coordinates. A summary of our conventions is 
given in the Appendix. In particular we use the follow- 
ing: 



{7^,7^} 
[7 M ,7l 



26^ U, 



(26) 
(27) 



where the j v denote the Euclidean Dirac matrices, and 
a^ v are the spin matrices. Since we are working in the 
chiral limit the Dirac operator is given by eq. I|12|) and is 
by definition hermitian: 



(28) 



The covariant derivative defines the field strength tensor 
in the fundamental representation as 



[V M ,V„] = -iF^. 



(29) 



The functional determinant of the fermions can be writ- 
ten as 



det(z^) = v/det 



(30) 



which following a method originally introduced by 
Schwinger can be further expressed as 



det(»y) = exp f-i Sp^ 



"^e-WWn ] . (3i) 



Here Sp is the functional trace. For its contribution to 
the effective action we have to properly normalize it, i.e. 
subtract the free zero-gluon part: 

Iogdet(t f) n = -5 S p/°° 7 - • ( 32 ) 

The square of JJ can be decomposed further, 

W 2 = 7/.7*V M V„ = - ({ 7 p,7.}V M V, + [ 7/t ,7„] V^Vu) ■ 

Since in the second term on the l.h.s. the commutator is 
antisymmetric we can also antisymmetrize 
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With Eqs. (I26l27l29l> one finds 

J7 2 = V 2 U + a^F^ . 
Equation l|32|) hence becomes 



(33) 



logdet(i^) n = --S P / — (34) 
2 Jo s 

x {exp [s(V 2 l4 + v^")] ~~ ex P I s d 2 !^] } ■ 

The 1-loop action is UV divergent. This comes from 
the fact that the running coupling constant is divergent 
at the tree level. Since QCD is a renormalizable theory, 
the tree level divergence has to be canceled by a 1-loop 
divergence. In order to control the divergent behavior 
of (|34J) we regularize the determinants by introducing a 
Pauli-Villars cutoff M in momentum space. This means 
that we use the so-called "quadrupole formula" : 



det(z y) r 



exp 



/det(- JJ 2 ) det(- ^ 2 + M 2 ) 
det(- @ 2 ) det(- f 2 + M 2 ) 



(35) 



ds 

— Sp 

s 



1 -e 



sM' 



The functional trace Sp can be taken by inserting any 
complete basis. We choose the plane wave basis: 



Spe~ sK = Tr / d 4 xlim f 

J y^x J (2tt 



-ip-y „-sK ip-x 



(36) 



where Tr is the remaining matrix trace over color and 
Lorentz indices. One can now drag the latter plane- wave 
exponent though the differential operator K until it can- 
cels with the former. This results in the shift of the 
derivatives inside the differential operator and in the fol- 
lowing representation of the functional trace poj ]: 

Spe- sK = Tr J d A x J J^L e -sK(d a ^ 9a+Va)l (37) 

The 1 at the end is meant to emphasize that the shifted 
operator acts on unity, so that for example any term that 
has a d a in the exponent and is brought all the way to 
the right, will vanish. According to l|37(l we now have 



logdet(^W=-^E J^sf 



k— — oo 



xTr {exp [s(V4 + iuj k ) 2 l4, + s(V, + ipif^A + sa^F^] 
-exp [s{iuj k ) 2 U + s(ipi) 2 U] } • (38) 
Let us now define 

B = Wi+iw k l 2 , (39) 
then eq. becomes 



logdet(i)7) n , r = -- J d 3 x J 



" ^Tre-^ 



1 - e 



k— — oo ' 



d 3 p 



(40) 



'){- 



+exp [(sB 2 + sV 2 + 2ispiVi)U + sa^F^] } 



This result is independent of the gauge group. In the fol- 
lowing we will work with SU(2). In particular we choose 
the background A4 fields to be a) static and b) diagonal, 
i.e. 



n 

2 ' 



then 



B = -i0(sc)y + iuj k l 2 . 



(41) 



(42) 



Here T3 is the third of the three Pauli matrices: 



n 



1 

1 



T-2 



-i 

1 



1 



• r3 = -1 H 43 ) 



V. COVARIANT DERIVATIVE EXPANSION 

A. Zeroth order - The fermionic potential 

In this order we set the spatial covariant derivative to 
zero, i.e. V* = 0, but sum over all powers of V4: 



k— — oo 



44) 



xTre" 



This can be evaluated explicitly. With eq. I|42|) it is easy 
to check that 



Tre 



sB 2 u _ 



e -i S (0-2w fc ) 2 +e -I s (0+2^ fc ) 2 



(45) 



where the factor 4 comes from Trl4. Since the fermionic 
energies are given by uj k = (2k + 1)ttT we can rewrite the 
terms in eq. l|43|l as 



e -\s(2uj k ±4,) 1 = e -s[2TrTk-<l> ± ] 2 



where we defined 



* 2 



(46) 



(47) 



The summation over w k and the integrations over s and 
p can now be performed along the lines of the bosonic 
case ([13), using the formula 



E 

k— — o 



3 7 

^_ e -s(27rTk~4>) 2 -sp 2 . 



1 U Id ^ 

= —log ch— — cos— 
\ T T 1 



The result is the following: 



(48) 



[logdet(i 



1(0) _ 



1 



12tt 2 T 



d 3 3 



(49) 



x [0 2 + (2nT - <M 2 + 0-(27rT - 0_) 2 - 2(^T) 4 ] mod ^ 



mod 4ttT 
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FIG. 1: The gluon (solid) and fermion (dashed) potentials for 
N f = 1 and -1 < v < 1. 



This potential is symmetric around = and is peri- 
odic with period AirT in contrast to the gluon potential 
which has period 2nT. The curvature around 
4> = of eq. H49|) gives the fermionic contribution to 
—m 2 D /T. Indeed we find — (7V/T)/6 which is in accor- 
dance with the known 1-loop result for the Debye mass 



(50) 



If we add the corresponding result from the gluons (see 
e.g. 0) 

1 



12tt 2 T 



<Tx 6 (2nT- 



od 2ttT 



(51) 



then we get the full result for two colors, namely (2 + 
N f /2)T/3. Introducing the variable v = ^/{2-kT) we 
have 4>± = ttT(v ± 1) and eq. (|49l) becomes 



2 /> 

[logdet(iX7)^] (0) =T 3 ^- J [d 3 x(l 



,2\2 



^] mod 2 ' 

The potential is then given by the (in the chiral case 
identical) contributions of all Nf quark flavors: 



(2?r) 2 r, o x9 , 

-Nf ^-r- 1 - ^ 2 2 - 1 , 

J 24 L J mod 2 



(53) 



This result is of course well known and can for example 
be found in the Appendix D of Q or in [l^. It can be 
compared to the pure Yang-Mills potential (0)01) 



i/YM _ ( 27r ) ,.2 ft ,,\2\ 

v — — 5 — v (1 — V) | mod 



(54) 



Both potentials are shown in Fig. ^ 

We see clearly that the YM potential has N c minima, 
which are the Z(N C ) symmetric points of = 0. Since 
there is no center symmetry for the fermions, we do not 
find the same situation. Indeed the fermion potential 
has a minimum at v = and a maximum at v = 1. Its 
period is doubled relative to the YM potential. This fact 
comes solely from the fundamental representation of the 
fermions. 



B. Leading terms in electric sector 

For an effective theory we are interested in the leading 
terms in the electric sector. We would like to stress that 
we are keeping all powers of the background field in 
our approach, but make an expansion in the spatial co- 
variant derivative. This means that we allow for an arbi- 
trary amplitude of the A4 fields but we assume that all 
the background fields are slowly varying and have mo- 
menta p < T. For the leading terms we hence expand to 
quadratic order in V^. Just as for the gluons and ghosts 
[l8| the technique is to expand 

Trexp {s [(B 2 + V? + 2ip l V l ) 1 4 + o^F^\ } (55) 

in powers of Vi using the following two master formulas 
for two non-commuting matrices A and B: 



dae aA Be (1 - a)A (56) 

pi pi — a. 

da dPe aA Be l3A Be^- a -^A + ... 

Jo Jo 



and 



[B,e A ] = f 
Jo 



= I d 1 e' yA [B,A} 



,(1-7)A 



(57) 



Since in eq. (|56|l powers of B are brought down in the 
expansion we identify B with the combinations of covari- 
ant derivatives in eq. H55|) and A is the rest. The electric 
field is identified as 



[V;,fi] = [V,,V 4 



-iFu 



-iEi 



-iEf T a 



(58) 



To the second order in Vi there are three terms con- 
tributing: 



Ti = 

T 2 = 

To = 



S Tr / da e asB2u (V 2 1 4 ) e ^ aB '^ 
Jo 

„2„2 /•! /-1-a 



(59) 



-— Tr / da 

3 Jo Jo 



x e 



s 2 Tr 



x er 



(Vil 4 je 

1 fl—a 
da / JR c <xsB 2 U 

(i Jo 



d(3e asB l4 (V,l4)(60) 

(l-o;-/3)sB 2 l4 



d(3e c 



(cr a 0F a(3 ) (61) 



The terms T\ and T2 are of the same structure as in the 
gluon case and can be shown (see Appendix B of (lq) to 
yield two gauge invariant contributions: 



1 all — a) 

2 V ' 



h = s 3 da j-~ +a(l - a) + ^sp 

xTre (i-«) s B 2 {js,Ei} e asB2 {B,Et}U , 
h = s 2 (~ - ^sp 2 ^Tre s£?2 (2i; 2 + i {B, [V 4 ,^]}) 1 4 . (63) 



(62) 
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In contrast to the pure Yang-Mills calculation from [l^ 
the Lorentz structure yields a factor 4 from Trl 4 and 
since we are dealing with fermions all matrices are in 
the fundamental representation. The second term in I2 
contains an anticommutator of B and the covariant diver- 
gence of the electric field, which is zero if the background 
field obeys the equation of motion. We will discuss it sep- 
arately in the next section and leave it out for the time 
being. 

What has to be evaluated is 



Here we used that io k = (2fc + l)irT and = 2-kTv. The 
function <&{v) is given by 



$0) = 2 



l + v 



1 - V 



(72) 



Here ip is the digamma function 



^(*) = |-bgr(«). 

az 



(73) 



(2) 



(64) 



The parameter /i is the UV-cutoff in divergent series: 



k— — oc 

-N F (Li + L 2 + L 3 



d 3 p 
J2nf 



f] C 2 

-e~ s P (h+I 2 +T 3 ) 
s 



For the summation over the Matsubara frequencies it 
turns out to be necessary to define a region of definition 
for cp. In particular we shall rescale this field as = 2ttTv 
and look at the interval — 1 < v < 1. In different regions 
of <\> the results will have different functional forms. We 
already saw for the fermion potential that it is symmetric 
in the interval — 1 < v < 1, and outside this region one 
has to continue analytically. 

We will start with the term T3. From 



(65) 



and 



Tr7 a 7 ( 37 7 75 = 4(<5 Q/ 3<5 7( 5 - 5 ai 5ps + S a sSpj) (66) 

it follows that 

Tr <7 a f3a 1 5F a f3F 1 5 = 2 F af3 F a p . (67) 

Since F a pF a p contributes to the electric sector as 2E k E k , 
eq. lj6*T)) is equal to 

T 3 = 4 S 2 Tr Cda C~ °df3 e^E^E^ 1 -^^ 
Jo Jo 
pi 

= 2s 2 Tr / dae asB2 E k e {1 ~ a)sB2 E k , (68) 
Jo 

where we used in the last line that the integrand is sym- 
metric in a and (3, invariant under a — » (1 — a) as well 
as the cyclicity of the trace. After integration over a, p 
and s we find the following structure for L3: 



1 

4^2 



((Elf + (Eff) l Sl + (E 3 f^S 3 } , (69) 



where 

00 I 

k— — 00 

00 

^3= E 

k— — 00 

1 



(j)U) k 7TJ 



1 



1 



4ttT 



\<f)-2uj k \ \(j} + 2LU k \ 
[4 ( 7 e- log + 



(71) 



00 ^ 

E^-E^g 



fc=l fc=l 
and is related to the Pauli-Villars mass as 

M 



fi, 



4ttT 



(74) 



(75) 



This subtraction scale for the running coupling constant 
has been known previously [|| and was also obtained in 

Next we will turn to the invariant I\ . After integration 
over a, p and s we find the following structure for L\: 



1 

4^2 



({Elf + (Eff) ^-S 1 + (E?f^S 3 (v) 



where S3 is given by eq. i|71|) and 

W> 2 +4^ 2 ) 



Si 



E 

k— — oc 



4> 3 u k 

4> 2 - 2<\>uj k + 4ujI 4> 2 + 24>u k + AujI 



(76) 
(77) 



2uj k 



2uj k 



= — [2 (log/i - log 4 - 2 7E ) - $(«/)] . 

Finally we investigate /2- Again after integration over 
a, p and s, L2 is of the form: 



L? 



1 

"4^2 



((i? 1 ) 2 + (Eff + (£?) 2 ) JS a 



(78) 



with 53 given by eq. (|T1|> . 

Collecting all terms from £1,2,3 we find the following 
results for the kinetic energy in the electric sector 



(2) 



[Sf-loopJ E 



d 3 x 
d 3 x 



((Elf + (E?f)fl F \v) + (EtffP(v) 
(EfA%f , ¥) 



f d x F-F-f^hA I WW 
J—yF l E l } 1 (v) + / 2 



, (79) 



where (v) = /a (v) — /{ (v) ■ In the second line of 
eq. lf7§)l we used 



(Eff - < W 



^14^-4 



(80) 
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cutoff M by A and find 



N 



S 



F. 



F, 



(P), 



3 (") 



24tt 2 



log 



A 



o7E 



4ttT 



(87) 



96tH 



-41og ■ 



A 



N 



f 



96tt 2 

with $(^) given by eq. i|72|) 



4ttT 

[4 ( 7E + log4) + $(^)] 



(89) 



FIG. 2: The functions f[ F) (solid), /^ F) (dotted) and / 3 (P) 
(dashed) without the UV divergent term. 



which follows from eq. (|41|l . The functions are given by 



A (F V) 



24tt 2 



log4^i, 



96tt 
N 

Wit 



2 [4 (te - log /i) + 
J [4 ( 7E + log4) + 



(81) 
(82) 
(83) 



We would like to stress once more that these functions are 
the results for the interval — 1 < v < 1. They are plotted 
in FigEI an d one sees that they are symmetric. To get 
outside the interval — 1 < v < 1 one has to continue 
analytically, and the functional form of the ft changes. 

We notice that the function f-f^ is constant, i.e. in- 
dependent of A4. However, it contains the UV divergent 
logfi that is necessary to renormalize the running cou- 
pling constant from the tree level action: 



F 2 



Fl 



ig 2 (M) 8 



—log 



M 

T 



11 

12 J 



7 Nr. 



-N 



(84) 



Here N c denotes the number of colors and Nf the number 
of flavors. We correctly obtained t he g luonic contribution 
to the charge renormalization in 18] . For the fermions 
the tree level divergence in the electric sector is 



EkEk 

WW) 



Nf M 
^24> g A- 



(85) 



From our result eq. (|81|l we find the correct UV divergent 
contribution 



EkEk 



Nt 

-E k E k — Mog 



M 



,7E 



(86) 



If we add the tree-level and the 1-loop action then the 
result should be UV finite. This is obtained by choosing 
the scale M in eq. I|86|) to be equal to the Pauli-Villars 
mass, which corresponds to the evaluation of the running 
coupling constant at the scale 47rT/exp(7E). In the ef- 
fective action we then have to replace the Pauli-Villars 



C. The "equation of motion" term 

In the previous section we left away the contribution 
of the second term in eq. I|t)3l) . Its contribution to the 
effective action is 



q N f f 3 r d z p r°° ds 



-sp 



(90) 



x 4s^ 



9 



sp' Tre sB (i{B,\yi,Ei]}) 



After integration and summation this becomes 



5] 



EM 



12tt 
Nf_ 
12tt 



, 3 TrQV;,^]^) 



(Fx 



ttT 
2ttT 



(91) 



Here we wrote the result once in terms of a covariant 
derivative in the fundamental representation and once 
through a covariant derivative in the adjoint representa- 
tion 

Df = diS ab + ! ach A\ and (A#i)° = Df e Ef, (92) 

in order to compare to the gluon results we 
obtained two terms: one comes solely from the non-zero 
Matsubara modes, while the other is the contribution of 
the zero mode alone. Our result here, eq. (|91|l . is equal 
to —Nf/2 times the first term of the gluon results. 

Eq. (|91fl is zero if the background field obeys the equa- 
tion of motion, D; L E; L — 0. Otherwise it depends on the 
behavior of A4 and Ei at spatial infinity. One can inte- 
grate eq. (|9"T)l by parts and gets 



N f 



'EM 



24tt 2 T 



^xiEfE^-diiEfAl)}. (93) 



It yields a contribution to the function eq. 1)87(1 plus a 
full derivative term. There are certain background fields, 
BPS dyons (23 being an example, where A4 — > const, 
and Ei ~ at spatial infinity. Therefore the full deriva- 
tive term in eq. I|93|l is nonzero. However, in the partic- 
ular case of the dyon, it satisfies the equation of motion 
and the two terms in eq. I|93|l cancel out. 
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D. Comparison to previous work 

In a related publication by Wirstam |28| an effective 
theory for QCD at high temperatures was derived by 
calculating gluon by gluon scattering at low momenta in 
terms of Feynman diagrams. In order to compare to the 
results of |23| we have to expand our functions I|87I89|) 
to quadratic order in v. For eq. (JH3 this gives naturally 
zero, and the remaining contribution is 



, 7C(3)iV/ 
' 384tt 4 T 3 



(94) 



which agrees with the result found in |28j if the gauge 
group is chosen to be SU(2). 



E. Leading terms in magnetic sector 

For an effective action in terms of magnetic fields we 
have to expand eq. (|40|l to quartic order in Vi . The basic 
idea of the calculation is, again, to use master equations 
(|56I58[) to drag covariant derivatives Vi to the right. One 
has for the commutators 



V*e 



is I dSe 
>o 



SsB' 



{B, Ei}e 



(l-S)sB 2 



V,V, e sB = < sB ' 



where 



[V.Vj ■, B 2 } = -zV,{6, E 3 } - t{B, Ei}Vj 



In this way one ultimately obtains gauge-invariant com- 
binations of the electric field in the fourth power, 
mixed terms containing both electric and magnetic fields, 
derivatives of the electric field and, finally, magnetic field 
squared. In this paper we restrict ourselves to the latter 
terms quadratic in the magnetic field Bi defined as 



D 1. ip r>a rpa 



(95) 



where Fj k — z[Vj,Vfc] in the fundamental representa- 
tion. For that reason we shall disregard the commuta- 
tors [Vi, A4] as they introduce powers of E^. In addition, 
we restrict ourselves to the magnetic field parallel to A4, 
i.e. Bi — B 3 T 3 . It means that we set the commuta- 
tor [Fy, At] = i (jS/i,Ej] — [S/i,Ej]) to zero. In practi- 
cal terms this means that we can drag all powers of the 
covariant derivative Vi as well as of the field strengths 
Fij through the exponentials of A4, as if they commute. 
Looking at the argument of the exponent in eq. H4U|) we 
see that terms which are quadratic in B, either do not 
contain the field strength tensor F^ v at all or consist only 
of powers of the latter. Mixing terms vanish upon inte- 
gration over momentum. Hence similar to the gluon case 



(see Appendix C of [ji|) we have to evaluate the follow- 
ing: 

[Sf-loop] m = ~ N f P°Sdet(< W)n,r] ( M ] (^) 

= —Nf [JVi - N 2 ] , 
where V\ only contains powers of the derivatives 



(2isY 
4! 



PiPiPfcPiVjVjVfcV; (97) 



(2is) 2 s 
3! 



PiPj [V 2 V 4 Vj + ViV 2 Vj + ViVjV 2 ]^ 



and V 2 comes from the field strength tensor alone: 



V 2 = s z Tr da d/3e ast> a ii F ii eP" s (98) 
Jo Jo 

x m Fi e {1 ~ a - p)sB2 

A V km r km~ 

We start with the evaluation of V\ . For the momentum 
integration in we use 



d 3 p 



d 3 p 



(2tt) : 



■PiPj 



d 3 p 

W) 



jPiPjPk p m e' 



-sp 



-sp 



-sp 



1 



1 1 

2^ (4tt S ) 3 / 2 UlJ ' 

1 1 
(2s) 2 (4tts) 3 / 



(99) 

(100) 

(101) 



X [SijSkm + SikSjm + Si m Sjk] , 

and obtain the following contribution to eq. 1961) : 

1 °° POO J 1 

^ = i^E I ^Tre^-lV^V,]^,^-]. (102) 

Since [V 2 ,V 3 ] 2 = -F^F^ = -2B k B k , where B k = 
B 3 . T 3 , this is equal to 



00 roo 



1 — ' r°° rt 

k— — 00 



and after integration over s and the summation over the 
Matsubara frequencies it becomes 

N ^-^JdMBf) 2 ^S 3 , (104) 

where S3 is again given by eq. H71|) . For V 2 we find, using 
eas. (|66l67[) : 

V 2 = 2 S 2 Trl da ( l ~dPe asB2 F l0 e? sB2 F.^ 1 -*-^ 2 . (105) 
Jo Jo 
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Since the contribution to the magnetic sector of FijFij is 
2B k B k this is equal to 

V 2 = ^Tifda [^"dpe^ 2 B k e 0sB2 B k e^- a - > B2 
Jo Jo 

= 2s 2 Tr (e sB2 B k B k ^ , (106) 

where we used the fact that the integrand is symmetric in 
a and /3, the cyclic property of the trace and eventually 
dragged the magnetic field to the right. One hence sees 
that V2 is of the same structure as V\ . Explicitly we find 
after the integrations over p, s and the summation over 
cuk that 



N 2 = 



±-Jd 3 X (B^S 3 , 



(107) 



where S3 is given by eq. . Adding eqs.JUHI HoTI) 

we find the following result for the kinetic energy in the 
magnetic sector: 



where the coefficient is given by: 



M F V) 



N 



f 



[4 (tb - log ji) 



(108) 



(109) 



96?r 2 

with $(1/) as in eq. 1(72)1, 

The function h\ above is the result for —1 < v < 
1 and it is symmetric in this interval. It also contains 
the necessary UV divergent contribution to cancel the 
tree-level divergence of the running coupling constant, 
eq. JEU: 



B k B k B k Bk. 
-log p = -— — r log 



24tt 2 



24tt 2 



M 



(110) 



Adding up the tree-level and 1-loop terms is obtained by 
replacing fx in eq. (|109|l by A/ (47rT)exp(7E), which corre- 
sponds to an evaluation of the running coupling constant 
at the scale 47rT/exp(7E). The final result is then: 



N 



f 



96tt 2 T 



-4 log 



A 
4^T 



•$(!/) 



(Ill) 



VI. CONCLUSIONS 

We have calculated the 1-loop contribution of massless 
quarks to the effective action at high temperatures for 
any value of A4 and hence of the Polyakov line. While 
we have a general result in the color-electric sector, we 
restrict ourselves to a magnetic field parallel to A4. The 
covariant derivative expansion of this action has the form: 



err J 



d 3 x 



-TV 



J5?ff>(„) 



(112) 



/ A a A a II 

where v = 2 *t 4 an< ^ i IS ^ ne ma g ne ti c field parallel 
in color space to A 4 . Because of the Bianchi identity, 
[-Fy,Ai] = i ([Vi,Ej] — [Vi,J5j]), in the case where the 
magnetic field is not parallel to A4 one also has to in- 
clude terms with electric field and its derivatives into the 
effective action, otherwise it will not be complete. 

The potential V F has double the period as compared 
to the gluon induced potential, is symmetric in v 
between -1 and 1 and has been known before. It has its 
minimum at \A&\ = and a maximum at \A^\ = ±27rT. 

The functions F$ and h[ F) given by eqs.JHTJ EHJ> and 
eq. (jlllf! are new. All functions, both in the electric 
and in the magnetic sector, are symmetric in v between 
-1 and 1, which reflects the fact, that fermions are 
in the fundamental representation of the color group. 
Our results can be used for studies of QCD at high 
but not infinite temperatures, where the Polyakov line 
experiences fluctuations which are large in amplitude 
but long ranged and where the dimensional reduction 
approach is too crude. 

Acknowledgments 

We are grateful to Chris Korthals-Altes for a critical 
reading of the manuscript and for helpful comments. 

Note added in proof 

After the submission, the paper by E. Megias, E. Ruiz 
Arriola and L.L. Salcedo, hep-ph/0312133 appeared on 
the net, in which the authors compute the coefficients 
in the parallel electric and magnetic sector. Apart from 
a constant our functions h\ and agree with their 
results. 



APPENDIX A: EUCLIDEAN COORDINATES 

In finite temperature QCD one needs the Euclidean 
formulation of path integrals in order to give the par- 
tition function the statistical-mechanics interpretation. 
Superscripts M will denote Minkowski coordinates, while 
superscripts E will refer to Euclidean coordinates. Note 
that throughout the paper we have used Euclidean coor- 
dinates without any explicit superscripts. 

For space-time coordinates we have 



,M 



,M 



For the gluon and fermion fields we have 



K = 



iA 



A; A 



M 



iip f 



The Dirac gamma matrices are related as: 



IX = 7cf 7f = ~*7f 7^ = 7 5 " 
The Euclidean gamma matrices are hermitian. 



(Al) 

(A2) 
(A3) 

(A4) 
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